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Figure SI Comparing statistics to detect and quantify introgression (A) window = 5 kb, 4Nr = 0.01 
Plots show means and standard deviations for D and f estimators, calculated over 100 simulated sequences 
(See Methods for details). Simulations covered 11 different values of f, the proportion of introgression. Gene 
flow was simulated either from P 3 to P 2 (conventional model, left-hand column) or from P 2 to P 3 (right-hand 



column). Dashed diagonal lines show the expectation of a perfect estimator of f. Two times of gene flow are 
shown: 0.1 (red) and 0.5 (blue) x 4N generations ago. 
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Figure SI Comparing statistics to detect and quantify introgression (B) window = 5 kb, 4Nr = 0.001 




Figure SI Comparing statistics to detect and quantify introgression (C) window = 5 kb, 4Nr = 0.1 
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Figure SI Comparing statistics to detect and quantify introgression (D) window = 1 kb, 4Nr = 0.01 
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Figure SI Comparing statistics to detect and quantify introgression (E) window = 1 kb, 4Nr = 0.001 




Figure SI Comparing statistics to detect and quantify introgression (F) window = 1 kb, 4Nr = 0.1 
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Figure SI Comparing statistics to detect and quantify introgression (G) window = 10 kb, 4Nr = 0.01 
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Figure SI Comparing statistics to detect and quantify introgression (H) window = 10 kb, 4Nr = 0.001 




Figure SI Comparing statistics to detect and quantify introgression (I) window = 10 kb, 4Nr = 0.1 
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Figure S2. Effects of recombination rate one the D statistic and f estimators. Boxplots show the 
distribution of values for 100 simulated 5 kb windows, with f = 0.1, three different values for the population 
recombination rate (4Nr) were tested: 0.001, 0.01 and 0.1. 
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Figure S3. Effects of genetic diversity on the D and f statistics in Heliconius whole genome data. 
A,B. Values of D and the three f statistics for non-overlapping 5 kb windows across the genome, plotted 
against nucleotide diversity, f values are only plotted for windows with D>0. Data from Martin et al. 2013. 
Taxa used are as follows, Pi: Heliconius melpomene aglaope, P 2 : Heliconius melpomene amaryllis, P3: 
Heliconius timareta thelxinoe, O: four Heliconius species from the silvaniform clade. Colored points show 
windows located within the wing patterning loci HmB (red) and HmYb (yellow), see Methods. 
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Figure S4. Variance in D and f statistics plotted in bins of nucleotide diversity 

Each statistic was binned into one of four bins according to average nucleotide diversity (tt) in the three 
populations. Variance was calculated for each bin. For all four statistics, only windows with D > 0 were 
considered because f statistics are meaningless when D is negative, and can range from 0 to negative infinity. 
For f G and fhom, windows with values greater than one were excluded for calculating variance. Four window 
sizes were examined: 5 kb (A), 10 kb (B), 20 kb (C) and 50 kb (D). 



Table SI. d X y between P 2 and P 3 for different sets of windows, for models with GENE FLOW FROM P 3 
to P 2 . 

The first three columns give the split times for the model and the time of gene flow in units of 4N 
generations ago (see Fig. 4A). The 'Alt' and 'Bg' columns give dxv between P 2 and P 3 for the 1000 Alternate 
and 9000 Background windows, respectively. For each of the introgression statistics, there are three columns 
giving d X y between P 2 and P 3 for outliers (windows with the top 10% of values), non-outliers (the remaining 
90% of windows) and then the percentage of the outliers that corresponded to real Alternate windows in the 
simulation (i.e. the level of accuracy of the statistic in identifying gene flow windows). Highlighted blocks 
indicate cases where d X y in the Alternate or outlier windows was significantly lower than in the Background 
or non-outlier windows. 
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Table S2. d X y between P 2 and P 3 for different sets of windows, for models with GENE FLOW FROM P 2 
to P 3 . 

The first three columns give the split times for the model and the time of gene flow in units of 4N 
generations ago (see Fig. 4B). The 'Alt' and 'Bg' columns give d XY between P 2 and P 3 for the 1000 Alternate 
and 9000 Background windows, respectively. For each of the introgression statistics, there are three columns 
giving d X y between P 2 and P 3 for outliers (windows with the top 10% of values), non-outliers (the remaining 
90% of windows) and then the percentage of the outliers that corresponded to real Alternate windows in the 
simulation (i.e. the level of accuracy of the statistic in identifying gene flow windows). Highlighted blocks 
indicate cases where d X y in the Alternate or outlier windows was significantly lower than in the Background 
or non-outlier windows. 
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Table S3. d X y between P 2 and P 3 for different sets of windows, for models with ANCESTRAL 
STRUCTURE 

The first three columns give the split times for the model in units of 4N generations ago (see Fig. 4C). The 
'Alt' and 'Bg' columns give dxv between P 2 and P 3 for the 1000 Alternate and 9000 Background windows, 
respectively. For each of the introgression statistics, there are three columns giving d X y between P 2 and P 3 for 
outliers (windows with the top 10% of values), non-outliers (the remaining 90% of windows) and then the 
percentage of the outliers that corresponded to real Alternate windows in the simulation (i.e. the level of 
accuracy of the statistic in identifying structure windows). Highlighted blocks indicate cases where dxv in the 
Alternate or outlier windows was significantly lower than in the Background or non-outlier windows. 



t?3 tl2 tsTR 


Alt Bg 


D 


fd 


fa 


fhom 


non- % 
outlier . |. 

outlier cor. 


non- % 
outlier u . 

outlier cor. 


non- % 
outlier . |. 

outlier cor. 


non- % 

outlier 

outlier cor. 


0.4 0.2 0.6 
0.4 0.2 0.8 
0.4 0.2 1 


0.018 0.018 
0.018 0.018 
0.018 0.018 


0.017 0.018 52.8 
0.017 0.018 80.2 
0.017 0.018 92.4 


0.017 0.018 58.9 
0.017 0.018 82.9 
0.017 0.018 93.5 


0.017 0.018 61 
0.017 0.018 82.7 
0.017 0.018 92.3 


0.017 0.018 63.4 
0.017 0.018 85.6 
0.017 0.018 94.7 


0.4 0.2 1.2 


0.018 0.018 


0.018 0.018 97 


0.018 0.018 97.9 


0.017 0.018 95.5 


0.018 0.018 97.9 


0.4 0.2 1.4 


0.018 0.018 


0.018 0.018 98.9 


0.018 0.018 99.2 


0.018 0.018 97.2 


0.018 0.018 99.3 


0.4 0.2 1.6 


0.018 0.018 


0.018 0.018 99.5 


0.018 0.018 99.7 


0.018 0.018 98.4 


0.018 0.018 99.6 


0.4 0.2 1.8 


0.018 0.018 


0.018 0.018 99.8 


0.018 0.018 100 


0.018 0.018 98.9 


0.018 0.018 99.8 


0.4 0.2 2 


0.018 0.018 


0.018 0.018 99.8 


0.018 0.018 99.9 


0.018 0.018 99 


0.018 0.018 99.8 


0.6 0.2 0.8 
0.6 0.2 1 


0.021 0.022 
0.022 0.022 


0.021 0.022 47.2 
0.021 0.022 73.2 


0.021 0.022 65.6 
0.021 0.022 86.7 


0.021 0.022 70.9 
0.021 0.022 90 


0.021 0.022 71.2 
0.021 0.022 90.9 


0.6 0.2 1.2 


0.022 0.022 


0.021 0.022 87.8 


0.021 0.022 96 


0.021 0.022 96.9 


0.021 0.022 97.2 


0.6 0.2 1.4 


0.021 0.022 


0.021 0.022 93.6 


0.021 0.022 98.5 


0.021 0.022 98.6 


0.021 0.022 99.1 


0.6 0.2 1.6 


0.022 0.022 


0.021 0.022 96.6 


0.022 0.022 99.4 


0.022 0.022 99.6 


0.022 0.022 99.6 


0.6 0.2 1.8 


0.022 0.022 


0.021 0.022 98 


0.022 0.022 99.7 


0.022 0.022 99.8 


0.022 0.022 99.9 


0.6 0.2 2 


0.022 0.022 


0.022 0.022 98.6 


0.022 0.022 99.9 


0.022 0.022 100 


0.022 0.022 100 


0.6 0.4 0.8 
0.6 0.4 1 
0.6 0.4 1.2 
0.6 0.4 1.4 


0.021 0.022 
0.022 0.022 
0.022 0.022 
0.021 0.022 


0.020 0.022 44.1 
0.021 0.022 70.8 
0.021 0.022 86.1 
0.021 0.022 92.5 


0.020 0.022 55.8 
0.021 0.022 80.4 
0.021 0.022 92.3 
0.021 0.022 96.9 


0.020 0.022 59.3 
0.021 0.022 80.9 
0.021 0.022 92.7 
0.021 0.022 96.6 


0.020 0.022 59.7 
0.021 0.022 82 
0.021 0.022 93.3 
0.021 0.022 97.3 


0.6 0.4 1.6 


0.022 0.022 


0.021 0.022 96.1 


0.022 0.022 98.8 


0.021 0.022 98.3 


0.022 0.022 98.8 


0.6 0.4 1.8 


0.022 0.022 


0.021 0.022 97.5 


0.021 0.022 99.4 


0.021 0.022 99.1 


0.021 0.022 99.4 
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0.025 0.026 72.6 
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0.025 0.026 88.1 


0.025 0.026 91.5 
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0.8 0.4 1.4 


0.025 0.025 


0.025 0.026 81.1 


0.025 0.025 95 


0.025 0.025 96.1 


0.025 0.025 96.4 


0.8 0.4 1.6 


0.025 0.025 


0.025 0.026 87.5 


0.025 0.025 98 


0.025 0.025 98.2 


0.025 0.025 98.4 


0.8 0.4 1.8 


0.026 0.025 


0.025 0.025 91.8 


0.026 0.025 99.5 


0.025 0.025 99.3 


0.026 0.025 99.6 


0.8 0.4 2 


0.025 0.025 


0.025 0.025 95.1 


0.025 0.025 99.8 


0.025 0.025 99.7 


0.025 0.025 99.7 


0.8 0.6 1 
0.8 0.6 1.2 


0.026 0.025 
0.026 0.025 


0.024 0.026 38.5 
0.025 0.026 64.9 
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0.025 0.026 96.8 
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0.025 0.026 97.1 


0.8 0.6 1.8 


0.026 0.025 


0.025 0.026 94 
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0.025 0.026 99 
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0.025 0.026 96.3 
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0.029 0.029 82.1 
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0.032 0.033 48.1 
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0.032 0.033 88 
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0.033 0.033 96.9 
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0.033 0.033 


0.032 0.033 53.3 


0.032 0.033 83.9 


0.032 0.033 86.1 


0.032 0.033 86.4 


1.2 0.8 1.8 


0.033 0.033 


0.032 0.033 69.7 


0.033 0.033 94.3 


0.033 0.033 94.8 


0.033 0.033 95 


1.2 0.8 2 


0.033 0.033 


0.032 0.033 81.2 


0.033 0.033 97.9 
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0.033 0.033 
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0.032 0.033 81.6 
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Table S4. d X y between P 2 and P 3 for different sets of windows, for NULL MODELS 

The first two columns give the split times for the model in units of 4N generations ago (see Fig. 4D). For 
each of the introgression statistics, there are two columns giving d X y between P 2 and P 3 for outliers (windows 
with the top 10% of values), and non-outliers (the remaining 90% of windows). Highlighted blocks indicate 
cases where d X y in the outlier windows was significantly lower than in the non-outlier windows. 
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